o 
o 

(N 



Multi-Particle Decays of Light Mesons Measured by 
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The PHENIX experiment at RHIC measured Kg 77, and tj-meson production at high pr 
in p+p, d+Au and Au+Au collisions at yijvjv = 200 GeV. Measurements performed in 
different hadronic decay channels give consistent results. This paper presents measured 
meson-to-7r° ratios and Nuclear Modification factors in the most central ci+ An and 
Au+Au collisions. No suppression seen in d+Au interactions is in contrast to a strong 
suppression of meson yields revealed in central Au+Au collisions at the same energy. 
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1 Analysis 

The layout of the PHENIX detector [2] and the decay modes of the particles presented in this 
analysis are shown in Fig. [T] The reconstruction begins with pairing photons (straight lines) 

and selecting the 7r° candidates 
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based on the invariant mass of 
the pair. The 7r° candidates are 
then combined between them- 
selves, other photons or with 
the charged tracks and corre- 
sponding invariant mass distri- 
butions are analyzed to extract 
the particle yields by simulta- 
neous fitting of the peaks and 
the background. The posi- 
tions of the mass peaks were 
found to be in agreement with 
the particle masses measured 
in vacuum and the widths of 
the peaks, depending predom- 
inantly on the detector resolu- 
tion, change from 10 MeV/c 



for 77 - 
for uj 



Figure 1: The PHENIX detector layout and the and the 
decays modes of studied particles. PHENIX subsystems 
not used in the analysis are shown gray. 



to 20 MeV/c^ 
and A'^ 



tt^tt" and to 30 MeV/c^ for lo 
tt'^j. The values above vary 
within less then 5 McV/c^ de- 
pending on the pt bin which 
agree with the widths resulting from the detector resolution. 

The analysis discussed in these paper is based on the event samples accumulated during 
PHENIX physics Runs3,4, and 5 with integral statistics, after quality assurance selection. 



*for the full list of authors see [3] 
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corresponding to the integrated luminosity of 1.5 nb^^ (p+p) 129 /ib^^((i+Au) and 2.5 pb^^ 
(Au+Au) in these runs respectively. The background conditions, depending on the mode of 
study for pT >4-5 GeV/c, is smaller than 1:5 to 1:20 in p+p and d+Au and 1:70 in Au+Au. 

The raw yields have to be corrected for the limited detector acceptance, the 7-trigger effi- 
ciency, various analysis cuts, the gamma conversions in the detector preceding the calorime- 
ter and the branching ratios of the specific decay mode. The magnitude of the corrections 
is calculated based on the full detector simulation and analysis of the data. The efficiencies 
measured in the detector configuration used during the p+p data taking are shown in Fig.O 
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Figure 2: Shown on the left is the geometrical acceptance of the PHENIX detector for various 
decay modes. The right panel shows the probability that a photon coming from meson decay 
fires PHENIX high pT 7-trigger. A unit corresponds to the standard "Minimum Biased" 
PHENIX trigger registering 23±2.2 mb. 

The phase space density distribution of the true three-body decay modes tt'^tt+tt^, known 
from the literature [6] was taken into account in the simulation. More details about this 
analysis can be found in [31 131 E] . For Au+Au events we applied an additional correction 
for reconstruction losses due to the detector occupancy. 

The systematic errors of the measurement are listed in Tab. [TJ The main source of error 
is the extraction of the raw yield made by fitting. The procedure is described in [3] . 



Source 


p + p d + Au 


UJ 7r'-'7 
p + p d + Au Au + Au 


p + p d + Au 


Acceptance 
Trigger efficiency 
Yield extraction 
MB trigger 


5-10 9-12 
3-10 5-7 
5-25 10-15 
10 8 


10- 20 8- 12 14- 16 

2-7 5 

5-15 10 15-35 
10 8 4 


10-25 10-20 
2-10 5 
7-30 9 
10 8 


Total 


15 - 25 18 - 22 


15 - 25 17 - 20 20 - 45 


20 - 40 15 - 25 



Table 1: Systematic errors (in %) for different decay modes and collision systems. Values 
with a range indicate minimum and maximum error in the px range of the measurement. 
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2 Results 



The results of the multi-particle decay 
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Figure 3: Invariant yields of K (top), rj (mid- 
dle), and ui (bottom) measured in p+p, p+A\i 
and Au-|-Au collisions at y/sNN=200 GcV. 
The solid line is the parameterized 7r° spec- 
tra measured in p -\- p [7] . The dashed lines 
arc the 7r° spectra scaled by the meson-to-7r° 
ratio in p + p and the number of the binary 
collisions. K+K~ data is taken from [9j 



measurements arc presented in Fig. [3] 
for K-mcson in the top panel, r/-mcson in 
the middle panel and w-meson in the bot- 
tom panel. 

The w-meson is measured in two decay 
modes, lo — s- tt'^tt^tt", uj 7r°7 and ?/- 
meson is measured in 77 — s- tt'^tt+tt", t] — > 
77. The results of the measurements for the 
same meson agree. The Kg tt^tt^ cannot 
be compared to Kg — > tt+tt" from PHENIX 
because of the detector-induced background 
at the mass of the K-meson. The results 
are in agreement with the STAR experiment 
measurement in tt+tt^ I^. Wc also see very 
good agreement between the results ob- 
tained in PHENIX physics Run3 and Run5 
and between the results measured in trig- 
gered and Minimum Biased event samples. 
In the latter case the correction shown in 
right panels of Fig. [2] docs not apply. 

In the p + p data the px range of the 
measurement is limited by the detector ac- 
ceptance on the lower side of the range and 
by the available statistics on the upper side. 
In Run5 the w-meson production in p + p is 
measured out to 13 GeV/c making it the 
second in p^-reach identified particle after 
the TT*^. In the heavier collision systems 
the combinatorial background effectively re- 
duces the available statistics. 

The lowest data points shown in each 
panel of Fig. [3] are the yields measured in 
p + p collision. Plotted above them are the 
Minimum Biased and 0-20% central c?-|-Au 
events data. The 60%-92% peripheral, Min- 
imum Bias and 0-20% central Au-|-Au col- 
lision spectra are plotted on top. For Kg 
such data are not available. All measure- 
ments are done at y/s, y^SNN = 200 GeV. 
Central d+Au and Au-t-Au data are scaled 
by different factors for the clarity. 

The solid line shown in each panel is the 
parameterization of the invariant yield of tt'' 
measured in p + p collisions [7] . For the r]- 
and the i^S-mesons this line is above the 
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data points and for the tj-mcson the points are much closer. Using this parameterization 
the non-identical meson ratios can be calculated. We find that in the p + p collisions these 
ratios are flat above the pT=2.5 GeV/c. Fitted by a constant the particle ratios are: w/tt^ 
= 0.81±0.02±0.07, rj/n"^ 0.48±0.02±0.02 and K°/tt° = 0.45±0.01±0.05. 

Each set of points shown in Fig. [3] has a corresponding dashed line. These lines are 
constructed in the following way: the 7r° spectra measured in p + p (solid line) is scaled with 
the corresponding meson-to-7r° ratio given above. Since all ratios are found to be flat in 
the region of the measurement, the scaled tt*^ reference corresponds to the invariant yield of 
the meson in the p + p collisions. Those yields are further scaled by the number of binary 
collisions for each centrality bin in d+An and Au+Au presented in the figure. 

As one can see for all analyzed mesons the d+Au data in Minimum Bias and 0-20% 
most central events are very close to the dashed lines. The ratio of the two is the Nuclear 
Modification Factor and for the analyzed mesons it was found to be flat within the errors of 
the measurement. For the Minimum Bias event sample the RdA^^ are above unit but agree 
with 1 within the errors of the measurement. 

In the peripheral Au+Au collisions the nuclear modification factor is not very different 
from 1 as the dashed line lays close to the points. This is not so in the Minimum Biased and 
0-20% central Au+Au events. For w-meson wc find the Raa to be 0.4±0.15. The r^-mcson 
production in central Au+Au events is suppressed by a factor of 5 compared to scaled p + p 
reference. 
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